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Abstract 

Using the nonlinear Schrodinger equation (Gross-Pitaevskii equation), the dynamics of a macro- 
scopic wave packet for Bose-Einstein condensates falling through double slits is analyzed. This 
problem is identified with a search for the fate of a soliton showing a head-on collision with a hard- 
walled obstacle of finite size. We explore the splitting of the wave packet and its reorganization to 
form an interference pattern. Particular attention is paid to the role of gravity {g) and repulsive 
nonlinearity (no) in the fringe pattern. The peak-to-peak distance in the fringe pattern and the 
number of interference peaks are found to be proportional to g"^^"^ and UQ^^^f^/^, respectively. We 
suggest a way of designing an experiment under controlled gravity and nonlinearity. 
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Recently a great number of theoretical and experimental effort has been devoted to Bose- 
Einstein condensates (BECs). The trapping techniques can be used to control the typical 
size of localized states, such as a Gaussian wave packet, and the technique of Feshbach 
resonance can allow both the strength and sign of the nonlinearity to be varied. 

BECs show a variety of interference phenomena that reflect the matter-wave nature. 
Among them, the BEC analog of Young's interference experiments is the most notewor- 
thy. As demonstrated by Andrews et al, a pair of BEC wave packets separately located 
in double potential wells begins to diffuse after the release of the trap and produces in- 
terference fringe o^i^i^ . Since this pioneering experiment, numerous studies on interference 
have been devoted to creating a stable double-well trap to control the relative phase of split 
condensates^ i^i^i'^i^i^i'^" , except for some novel experiments based on an atom-based Mach- 
Zehnder interferometer that exploits a superposition of internal states of BEGii*i2iH. 

On the contrary, however, no work has been devoted to the BEC analog of single electron 
(photon) dynamics through real double slits (DSs). The release of a double- well trap cannot 
precisely mimic the passage through DSs. The former treats the second-order interference 
between two wave packets (WPs) with initially independent phases, while the latter treats 
the first-order interference between two fractions of the initially identical WP. A single BEC 
wave packet falling through DSs under uniform gravity will split into distorted pieces. How- 
ever the fate of these split pieces is not clear because of the complicated effect of diffraction 
at the slit openings (windows). Furthermore, the effects of gravity and nonlinearity on 
double-slit interference have not been examined quantitatively up to now. 

Both trapped and developing BECs are described using the Gross-Pitaevskii equation 
(CPE) or the nonlinear Schrodinger equation (NSE)^'^'i2iiiii^. In this letter, by using NSE 
(CPE), the dynamics of a macroscopic WP falling through the DSs is analyzed. In contrast 
to the ordinary quantum mechanics, BEC has advantage in that it shows a continuous 
time evolution of WP dynamics without being affected by the subtle problem of quantum 
demolition measurement. We here explore how WP splitting occurs at the slit openings and 
how the split pieces reorganize to form the interference pattern. We pay special attention 
to the effect of gravity and nonlinearity on the fringe pattern. We also discuss the designing 
of the experiment. 

At low temperatures, A^^o atoms occupy the single quantum state and the many-body 
wave function is described by a single macroscopic wave function (p which satisfies NSE, 



2 



i.e., the time-dependent Schrodinger equation with the additional mean-field type nonlinear 
term. In the presence of a harmonic trap and gravity, GPE (NSE) becomes 

i^-Q^ = "^^^'^ + ^t™p(r)0 + U gravity {T^)<j) + Mo|0| V, (l) 



where Uq is the nonlinear interaction defined by Uq = ^'^^ ° A^q [in 2-dimensional (2D) space] 
or ^'^^"' Nq (in 3D space) with a, m and / being scattering length, atomic mass and the char- 
acteristic length (to be defined below), respectively. Vtrapij) and Ugravityij) are potentials 
for a harmonic trap and gravity, respectively. 

Practically, a BEG WP in a harmonic trap exhibits a cigar-like shape.^^i^ When the WP 
lengthens in the z direction to form a nearly cylindrical shape, cross-sectional dynamics 
of the WP, which is observed on the x-y plane bisecting the z axis of the WP, can be 
approximately described using a 2D GPE without losing any physical meaning. Thus we 
perform 2D analysis that requires less computational cost than the 3D version. We assume 
an initial circularly trapped WP that falls toward the positive y direction under uniform 
gravity (gravity constant g). Then we set = ^ + ^ and 

Vtrap{Tc) = ^muj'^{x'^ + y'^)] Ugravityir) = -mgy. (2) 



Using the characteristic length / = ^JTiJJjmj), one can make all variables dimensionless as 

t' = cut, x' = j,y' = f , 0' = (3) 

Then GPE is reduced to 

1/92 52 \ , 1 



where Uq = AnNoa/l and g' = g/{luj'^) are also dimensionless constants. Initially, we ig- 
nore the contribution from gravity and, using the confining potential, prepare a circularly 
symmetric Gaussian WP with its center of mass at the origin r = (0, 0), 



^ini = : -771 exp . , (5) 



1/4 ^-'^Jr I 1 



which minimizes the total energy E' = J (i^r'(||V0'p + K^apl^'P + y'I^'H. Here, the de 
Broglie wavelength of an initially stationary WP can be longer than the size of the WP 
itself. Thus small variation of the WP shape contributes little to the global phenomena. 



The Gaussian WP profile is a simple substitute for an exact solution of GPE as far as 
short-time dynamics of the WP is concerned. 

Our strategy is as follows: (1) release the confining potential and switch on the uniform 
gravity; (2) WP falls on DS; (3) WP splits as it passes through the DS; (4) the split pieces 
reorganize to form a fringe pattern. Dynamics of the WP is obeyed by GPE in eq. (@]) together 
with the Dirichlet boundary condition at DS: hard walls forming the DS are taken to be one 
dimensional along the x direction. We assume the interference phenomena appear in the 
short-time dynamics shorter than the decoherence time and use GPE even after releasing 
the trap, closely following the procedures in most of the existing worka^i^jiSiiiii^. 

Given any arbitrary value -Uq, we systematically investigate the effect of gravity on the 
fringe pattern, starting from a common WP with fixed width. For this purpose, we always 
choose WP in the case of vanishing nonlinearity (mq = 0), 
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as an initial state for any value of nonlinearity. This means that we have replaced Vtrap in 
eq.(ID) by the Wg-dependent confining potential, + but have taken the scaling 

in eq.© using I proper for Vtrap with Uq = 0. Since we consider the falling of WP after 
releasing Vtrap, the formal change of the confining potential has no effect on the dynamics 
under gravity. 

We numerically integrate GPE in three ways, i.e., alternating direction implicit method 
(ADI method), Crank-Nicholson's method, and split-step method, all of which are found 
to yield identical results. Figure ^ shows the x'-y' plane of size 40x40 that we employ. 
Units for spatial and temporal grids are 6x' = 6y' = 0.1, 6t = 0.002. The center of WP, 
being located at the origin at t = 0, begins to fall along the y' direction due to the uniform 
gravity. DS is set at y' = 5. Lengths for the central stopper and a pair of windows are 
3 and 1, respectively. The fringe pattern will be observed on a virtual screen located at 
y' = Lsc = 18. The final time of the present computation is the time when the classical 



particle (center of mass) reaches the screen, and is given hy t = Tmax = y'^Lgc/g'- Up to 
this time, we detect little fraction of WP escaping outside the whole grid area through the 
boundaries at x' = ±20. 

The time evolution of WP in the linear case (uq = 0) is given in Fig. ^ The gravity 
is fixed a.t g' = 1 . Figure ^ shows that the falling WP breaks into a pair of pieces during 
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FIG. 1: Wave packet dynamics with u'q = 0, g' = 1. Intensity distributions of WP at t' = 0.0 and 
t' =1.0-6.0. Topmost panel illustrates the setting for double-slits experiments. 

passage through the DS and then begins to form a clear fringe pattern. Above the DS 
we find a remnant of WP, which has failed to fall through the DS, interferes parallel to 
DS. Because of the gravity, the remnant will also fall through DS sooner or later in due 
course. Our additional calculation (to be described elsewhere) indicates that corresponding 
WP dynamics in the case when either one of the DS windows is closed leads to no fringe 
pattern and merely to a single humped diffusive WP. To see the fringe pattern, therefore, 
it is essential to open both windows. The results are consistent with the result of a single 
electron interference using a DS. 

Taking the result for the linear case and g' = 1 a.s a standard pattern, let us proceed to 
the investigation of the role of tunable nonlinearity and gravity. We shall focus on the fringe 
pattern on the screen at y' = Lgc at t = T^ax- We choose u'q = —3, 0, 20 for the nonlinear 
term and g' =1, 3, 5 for gravity. Figure El which includes nine fringe patterns, represents 
the intensity distribution of WP at t = Tmax, and Fig. El shows their cross sections on 
the screen, with the total probability amplitude for this section being normalized to unity. 
We find that the nonlinearity is responsible for the localization or delocalization of WP, 
while gravity governs the velocity of WP arriving at the screen and is responsible for the 
de Broglie wavelength. For a fixed nonlinearity, the increased gravity will shorten the de 
Broglie wavelength at time t = Tmax- On the other hand, for a fixed gravity, the increased 
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FIG. 2: Intensity distribution of WP. From top to bottom: 151' =1, 3, and 5. From left to right: 



u'q = -3, 0, and 20. 
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FIG. 3: Fringe pattern on screen. Probability (0 — 0.2) versus horizontal coordinate. Arrangement 
of panels is the same as in Fig. [21 

nonlinearity u'q will broaden the fringe pattern. 

Let us develop a theoretical insight in greater detail. Under a fixed nonlinearity, by using 
Fig. El one can evaluate the peak-to-peak interval. The gravity dependence of the main 
interval (between the central and second peaks) is shown in Fig. |3 For each value of Mq, the 
interval decreases monotonically as a function of g' . The peak interval is proportional to the 
de Broglie wavelength defined by A' = In j-p' . If we choose p' to be the momentum for the 
classical particle corresponding to the center of WP, = Vmax = g'Tmax = g'\/'^Lsc/g' = 
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FIG. 4: Main interval versus gravity 



\/2g'Lgc- Then we obtain 



A' 



27r 



(7) 



regardless of the value u'q, the validity of which is confirmed in Fig. |31 

On the other hand, the nonlinearity is responsible for WP broadening, namely, the exten- 
sion of the area in which a fringe pattern can be formed. We shall make a simple evaluation 
of WP width at t = T„iax- Let us approximate the falling WP by a cylinder (0 = 7c~^^'^R~^) 
with radius R. By suppressing gravity, GPE becomes circularly symmetric and is given by 

1 d 



(8) 



2dR^ 2R OR 

Using the above cylindrical solution in this equation, we have R = i{uQ — TT/2)7i^^R~^, whose 



solution is i? ~ a/ {u'q — 7r/2)t in the large t region. Therefore, at t = T„ 

R.r oc 



'u'o-n/2 



(9) 

We then proceed to evaluate the peak number. The peak-to-peak interval Ax^^^i^ is 
approximately equal to the de Broglie wavelength A' at t = T^ax, so Ax^g^^, oc A' oc 1/y/g'. 
Eventually, the peak number (PN) at the screen can be given by 

R.,r If. vr 



PN 



Ax' 



oc 



peak 



Un - 
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(10) 



This approximate analytical expression for PN explains how the peak number depends 
on the nonlinearity and gravity. Table Q] obtained from our numerical fringe pattern on the 
screen, lists PN for various values of g' and u'q and is consistent with eq. p()|l . 

The transmission probability of WP at t = T^ax is the portion of WP that lies below 
the DS at that time. This probability is affected by the nonlinearity. Let us define Pr(y4) 
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TABLE I: Peak number 



and Pr(i?) for each portion lying above and below DS, respectively. They are calculated at 



f = T 



as 



Pt{A or B) 



dx'dy 



(11) 
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Noting Pr(y4) + Pr(i?) = 1, Pr(i?) gives the transmission probability. Pr(i?) at t = Tmax is 
found to depend on u'^^g' and takes values in the range 0.1 < Pr(_B) < 0.3. For larger Mq, 
WP broadens before reaching the DS and therefore Pr(i?) diminishes. On the other hand, 
for increased g\ WP reaches the DS without broadening very much because of the highly 
accelerated motion, and a large portion of WP dwells around the central stopper, which 
suppresses Pr(i?). Thus Pr(i?) is sensitive to the values of u'q and g' . However, qualitative 
features of the interference, such as g' and u'q dependences of the peak-to-peak distance and 
peak number, remain unchanged. 

We have so far investigated WP dynamics using scaled time and coordinates. Here we 
suggest an experimental situation to which the present result is applicable. Let us introduce 
the material-specific constant a = ^Jh/m. The value g' that we have employed determines 
uo and /. With the use of / 



^Jh/{muJ) and g' = g/iluj"^ 

2/3 / ^ ^ 2/3 



m g 
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and 



a 
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a J 



-1/3 



one obtains 

2/3 

9_ ' 
9' 



(12) 



(13) 



The last value represents the scale of, for example, the WP width, the DS, and the stopper. 
Let us consider more concretely the case of ^''Rb and g' = 1. Noting that a = 2.726 x 10^^ 
and the gravity constant on the earth g = 9.80665 m/s, we obtain / = 0.383 /im. In a 
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laboratory of the space station, one can expect a weak gravity, e.g., g = (?/1000. Then the 
experiment corresponding to g' = 1 is reahzed on a length scale / that satisfies 



from which we obtain / = 3.83 /im. The above equality indicates that on decreasing the 
gravity scale by three orders, the corresponding length scale increases by one order. In other 
words, the length scale is almost insensitive to the change of gravity and the size of DS can 
be kept almost constant when designing the present experiment under a weak gravitational 
field. 

We have explored the effect of nonlinearity and gravity on the macroscopic fringe pattern 
generated by BEC passing through DS. Most of the activities so far concerning the BEC 
analog of Young's interference experiment referred to a double-well trap leading to the 
2nd-order interference between two initially separated WPs with uncontrolled independent 
phases. By contrast we mimic the BEC analog of the Ist-order interference (a single electron 
interference). A study of the BEC dynamics through the DS is identified with a search for 
the fate of a 2D soliton showing a head-on collision with a hard-walled obstacle of finite 
size. In our simulation, we find that split pieces of WP below the DS successfully reorganize 
to form a clear fringe pattern. The analytic expressions are derived for the peak-to-peak 
distance and peak number as a function of gravity and nonlinearity. It should be noted that 
our preliminary simulation for a more realistic 3D version of the present system yielded the 
same result as found here. A further analysis of wave-function patterns, such as nodal lines 
and phase rigidity, which is developed in the context of quantum chaosi&, are subjects which 
we intend to study in the future. 
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